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ABSTRACT 


We  analyze  the  mechanism  by  which  a  particular  structure 
implodes  under  hydrostatic  pressure,  calculate  the  accelerated 
motions  of  the  boundary  surface,  and  calculate  the  magnitude  and 
distribution  of  the  sounds  which  are  produced.  We  then  describe 
experiments  to  verify  the  theoretical  analysis  of  the  motion 
and  the  sound  generation.  The  limited  and  incomplete  experi¬ 
mental  data  which  are  presently  available  substantially  confirm 
the  theoretical  analysis.  In  particular,  the  theory  and  experi¬ 
ments  show  that  the  principal  component  of  the  transient  sound 
is  the  same  in  all  directions  and  is  proportional  to  the 
second  time  derivative  of  the  volume  enclosed  by  the  wetted 
surface  of  the  imploding  structure.  There  are  some  uncertain¬ 
ties  in  the  interpretation  of  the  data,  but  these  should  be 
resolved  by  the  remaining  test  program. 

ADMINISTRATIVE  INFORMATION 

This  work  was  principally  performed  under  the  Foundation  Research 
Program  of  the  Naval  Ship  Research  and  Development  Center,  Task  Area 
ZR0110101,  Task  IR,  Work  Units  903-008  and  945-051.  Portions  of  the 
material  in  this  report  were  also  contained  in  a  thesis  submitted  by  W.T. 
Reader  to  the  Faculty  of  the  School  of  Engineering  and  Architecture, 
Catholic  University  of  America,  in  partial  fulfillment  of  the  requirements 
for  the  degree  Doctor  of  Philosophy  (May  1970). 
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INTRODUCTION 


Our  objective  la  (1)  to  explain  the  magnitude  and  distribution  of 
the  sound  which  la  generated  by  the  implosion  of  a  hollow  a true tu re  under 
external  preeaure  in  deep  water  and  (2)  to  verify  Che  theoretical  analysis 
by  experiments  under  controlled  conditions, 

By  an  implosion  wo  moan  a  sudden  collapse  of  the  structure  which 
occurs  when  the  external  pressure  reaches  soma  characteristic  critical 
value  and  triggers  an  unstable  deformation.  The  deformation  is  unstable 
because  the  external  pressure  forces  continue  to  exceed  the  internal 
restoring  forces  generated  by  the  deformation.  Eventually  the  forces 
balance,  the  kinetic  energy  is  dissipated,  and  the  motion  is  arrested, 
During  the  implosion,  sound  can  be  generated  by  a  variety  of 
processes  such  as  the  inrush  of  water,  the  Impact  of  two  parts  of  the 
structure,  or  the  sudden  compression  of  the  trapped  air.  These  processes 
all  generate  accelerated  motions  of  the  boundary  surface  of  the  structure 
and  thus  radiate  sound  into  the  water. 

respite  the  great  complexity  and  variability  in  possible  motions 
of  the  boundary  surface,  however,  we  show  that  for  a  wide  range  of  prac¬ 
tical  cases,  the  sound  pressure  in  the  far  field  of  a  uniform  medium  is 
the  same  for  all  directions  and  is  simply  proportional  to  the  second  timo 
derivative  of  the  volume  enclosed  by  the  wetted  surface.  Also,  we  show 
how  the  history  of  this  volume  acceleration  depends  on  the  Initial 
pressure  difference  and  on  the  internal  forces  which  resist  the  collapse. 

This  analysis  is  confined  to  predicting  the  shape  of  the  sound 
pulse  as  it  Is  affected  by  the  implosion  mechanism  and  as  it  would  be  in 
a  uniform  unbounded  medium.  In  any  real  situation  where  the  pulse  propa¬ 
gates  over  very  long  ranges  in  the  sea,  the  initial  shape  of  the  pulse 
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would  bo  enormously  distorted  by  multipatha,  by  velocity  gradients  in 
tho  wo* or,  ond  by  scattering!  the  pattern  of  implosion  aouud  at  aomo 
diatant  point  would  bo  determined  principally  by  the  geography  of  tho 
path, 

Tho  experiments  aro  mo do  with  a  simple  imploding  structure  for 
which  both  tho  coliopaing  motlona  and  tho  aound  pulao  eon  ho  predicted 
theoretically  and  moaaurod  with  proeiaion, 

Tho  prooont  roport  ia  only  a  progroaa  raport  bocauao  tho  expori- 
montal  program  haa  not  boon  oomplotod,  However,  ovon  tho  ilmitod  experi- 
montal  data  that  aro  available  eloarly  verify  tho  main  foaturoa  of  tho 
theoretical  analyaia,  Tho  additional  experimental  data  from  tho  remainder 
of  tho  program  aro  expected  to  rooolvo  any  uncortaintioa  remaining  in  tho 
onalyoia. 


THEORETICAL  ANALYSIS 
THEORY  AND  ANALYSIS  OF  THE  GENERAL  CASE 

Wo  atart  with  tho  Kirchhoff  equation1  for  tho  Inatantaneoua  aound 
proaauro  p(r,t)  at  an  arbitrary  field  point  r  and  time  t  in  torma  of  tho 
prior  valuoa  of  tho  aound  proaauro,  ita  time  derivative,  and  lea  normal 
gradient  at  all  pointa  i  on  n  fixod  aurfaco  which  encloses  the  imploding 
body. 


rr r |e 1 42. .  rrjii:*  <2.  rrM^sjw 
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Hare  n(a)  is  a  local  normal  pointing  outward,  R  -  r  -  j[,  and  the  bracket 
notation  signifies  that  the  term  within  the  brackets  is  to  be  evaluated 
^References  are  listed  on  page  43, 
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for  the  retarded  time  t  (^)  »  t  ■  P/o,  e,g,  , 


|p(|.0)  "  P<!>e  >  **  P<i*f  •  */«> 


(2) 


Thun,  the  entire  aound  preeaure  field  la  expressed  by  the  history  of  the 
aouud  preeaure  at  the  closed  aurface,  However,  Che  form  of  the  relation 
ie  inconvenient  to  analyte  (l)  becauae  the  retarded  time  at  which  each 
Integrand  la  to  be  evaluated  varlea  from  point  to  point  In  each  Integral, 
(2)  becauae  each  integrand  la  a  complicated  mixed  function  of  the  coor- 
dlnatea  of  the  field  point  £  aa  well  aa  the  aurface  point  a,  and  (3) 
becauae  the  point  j  identiflea  a  geometrical  location  on  a  fixed  aurface 
rather  chan  a  material  point  on  the  body  for  which  it  would  be  eimpler 
to  write  an  equation  of  motion, 

Hence,  we  take  the  origin  at  aome  convenient  reference  point  (to 
be  determined)  within  the  aurface  and  refer  all  terma  in  the  intagranda 
to  a  common  retarded  time  tQ  "  t  -  r/c.  We  do  this  by  expanding  each 

it  it 

function  of  the  retarded  time  aa  a  Taylor  aerlea  about  t  ■  t  ,  Thua  for 

o 

the  function  q(a,t  )  •  J^p/^n) 


q(«.t*) 


o)  dt*" 


★ 

t 


(3) 


and  similar  aerlea  expansions  apply  to  [p]  and  [p],  Also  we  restrict  the 
application  of  Equation  (1)  to  the  far  field  where  r  la  much  larger  than 
the  maximum  dimension  of  the  closed  aurface  o,  Then 


I  -  7  (1  +  0(e/r)] 


(4) 
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t*  -  t*  *  ■  *  «o«(jil)  U  +  0<«/r>  1  (*> 

Now  we  substitute  Squat  loot  (3),  <A> ,  and  ($>  into  Squat  ion  (1),  ret  ain 
only  those  tsrms  which  are  0(»/r),  and  apply  the  boundary  condition 


% 


-  pv 


(6) 


where  v(j[,t)  ia  th#  outward  normal  velocity  at  The  reault  la  a 

conventional  multipole  expansion  with  contrihutiona  from  the  firat  two 
integrate  in  Rquation  (l)  but  with  no  far*field  contribution  from  the 
third  integral, 


"<£•'>  ‘  (  X  II  [~ 


C0i(l.l)~|m  l  d"Vafrl 


J  m: 


dt 


m 


da 


<»> 


nv 

it 

where  all  terma  are  to  be  evaluated  at  the  conmon  time  t-  t  -  r/c, 

o 

We  tranaform  this  to  a  form  where  ail  Integrals  are  independent  of  the 
position  of  the  field  point  and,  furthermore,  where  each  multipole  term 
haa  a  simple  physical  interpretation. 


p<i.u  -  [  If  uW°] 

+  7i  [  II  (f*i  -  pa> ia]  ■  i 


+  TZ? 


£  •  ~-a  [  }  JJ  £<P*X  -  2p n)  da]  •  f  +  ... 


(8) 
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Tin*  bracketed  factor  tn  ths  first  tsrm  is  limply  pV  whir#  p$(t.*) 
is  tlu*  outward  flux  of  water  chrou|l\  ths  closed  surface,  Ths  aound  ft* Id 
(a  (ho  same  In  all  direct. Iona  and  la  tho  aama  aa  that  of  a  aimple  source 
(monopole)  whose  Instantaneous  source  strength  is  ^(tQ), 

The  second  term  Is  a  dipole  field  fr  cos(£,£)/Anrc  which  is  the 
same  aa  though  a  concentrated  force  £(tQ)  were  applied  at  the  origin  in 
an  unbounded  uniform  medium,  and  where  the  vector  force  £ 

The  second  term  in  £  la  simply  the  unsteady  force  exerted  by  the  water  on 
the  surface;  the  first  term  can  be  interpreted  as  an  inertia  force  asso¬ 
ciated  with  the  linear  momentum  of  the  surface  motion, 

The  third  term  in  Equation  (8)  is  a  quadrupole  field 
(l*Q‘.£)/(^nrcu )  which  is  the  same  as  though  an  equivalent  dyadic  torque 
•  j  II  ji (pvs  -  2pn)  do  were  applied  at  the  origin  in  an  unbounded 
medium.  The  second  part  of  Q  is  simply  the  unsteady  torque  exerted  by 
the  water  on  the  closed  surface;  the  first  part  is  an  Inertia  moment 
associated  with  the  angular  momentum  of  the  surface  motion. 

The  higher  multipole  terms  of  Eouation  (7)  can  also  be  written  in 
the  tensor  form  of  Equation  (8)  which  separates  the  direction  of  the 
field  point  from  the  surface  integrals, 

Tn  the  subsequent  analysis,  we  treat  £  as  the  Lagrangian  coordi¬ 
nate  of  a  material  point  of  the  wetted  surface  of  the  Imploding  body. 

This  does  not  invalidate  the  Kirchhoff  equation  provided  the  surface  does 
not  change  appreciably  during  the  time  interval  for  a  sound  wave  to 
transit  across  the  body,  i.e.,  provided  j  «  c.  Under  this  same  con¬ 
dition,  it  does  not  affect  the  Taylor  expansion  coefficient  given  by 
Equation  (5),  for 
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t*  -  t*  -  (t  -  R<t*)/c)  -  tc  -  r(t^>/c ) 


r<co>  "  R<0  R<0  "  R(t*> 

- c -  4 - c - 

a  £  C0*(liX>  +  (*  -  '*>  I  co“Ci*r) 

•’*  c*  "  Co  "  c  coi^l»^  [l  ‘  c  coaCk»l)] 

at  |  coa(s,r)  (9) 

In  practical  caaaa  of  implosion,  we  can  generally  expect  the  mono¬ 
pole  radiation  to  dominate;  our  experience  indicates  that  an  implosion  is 
accompanied  by  some  "crushing,"  i,e,,  some  volume  change,  and  Equation  (7) 
shows  that  a  monopole  aourco  is  the  most  efficient  multipole  sound  source 
because  of  the  c~m  factor,  However,  there  is  no  absolute  necessity  for 
the  higher  multipolo  terms  to  be  negligible,  and  indeed  it  is  simple  to 
invent  situations  whore  the  dominant  contribution  would  be  from  some 
higher  multipole.  But  the  external  pressures  that  initiate  the  implosion 
are  the  same  in  all  directions.  The  implosion  can  generate  sound  pressures 
with  a  higher  multipole  symmetry  only  if  the  internal  forces  and  the 

motions  of  the  wetted  surface  have  this  3ame  symmetry, 

THEORY  OF  A  SIMPLIFIED  IMPLOSION  DEVICE 

The  implosion  device  which  is  used  in  the  experiments  can  be 
idealized  for  analysis  purposes  in  the  simple  schematic  form  of  Figure  1, 

Initially,  the  piston  is  retained  at  the  top  of  the  air  chamber 
where  it  seals  off  the  internal  air  at  atmospheric  pressure  from  the 
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AltTON  IN  INITIAL  POSITION 


CAII 


SPNINO 


Flgurt  1  -  Idtaliisd  Implosion  Dovict 

external  water  at  some  higher  hydrostatic  pressure.  The  piston  is  then 
released  at  some  preselected  external  pressure  and  is  driven  down  the 
chamber  by  the  water  pressure  with  increasing  velocity  until  the  motion 
is  arrested  by  impact  with  the  spring  or  by  the  pressure  of  the  trapped 
air.  The  sharp  deceleration,  or  impact,  results  in  high  accelerations 
of  the  wetted  surface  and  in  a  significant  production  of  sound. 

To  analyze  this  motion,  we  assume  that  the  periods  of  any  elastic 
vibrations  of  the  piston  and  the  case  are  all  very  much  smaller  than  the 
time  intervals  in  which  the  motions  of  the  piston  and  case  change 
appreciably.  This  means  that  we  can  treat  the  piston  and  case  as  rigid 
bodies.  We  also  assume  that  the  transit  time  of  a  sound  wave  in  the  water 
across  the  case  is  very  much  smaller  than  the  time  interval  for  the  motion 
to  change  appreciably.  This  means  that  the  unsteady  pressure  distribution 
in  the  water  adjacent  to  the  piston  and  case  is  the  same  as  though  the 
water  were  incompressible,  i.e.,  the  same  as  though  the  water  acted  as 
an  added  inertia  whose  magnitude  depended  only  on  the  patterns  of  motion 
of  the  pi3ton  and  the  case. 


zzzzzzzzzz 
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Thus,  if  Ax  ■  and  Aa  -  naa  are  Che  section  areas  of  piston  and 
case,  respectively,  and  if  x:(t)  and  xB(t)  are  the  upward  displacements, 
then  the  kinetic  energy  of  the  water  can  be  written  in  the  form 

Lxft?  La*a 

T-— +  _  (10) 

where  the  added  water  mass  due  to  piston  motion  is 

1*1  •  2paf  +  oA^Xg  -  xx )  (11) 

The  first  term  is  the  added  mass  for  a  thin  disk  of  radius  a  when 
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one  face  moves  with  uniform  velocity  and  the  other  face  is  stationary. 

It  thus  approximates  the  added  mass  of  a  semi-infinite  cylinder  moving 
end-on.  The  second  term  is  the  mass  of  water  in  the  tube.  To  the  same 
approximation,  the  added  water  mass  due  to  the  motion  of  the  case  is 

La  *  2pa|  (12) 


These  assumptions  also  imply  that  the  upward  force  on  the  water  due  to 
the  motion  of  the  piston  is  (from  Lagrange's  Equation) 

(13) 

(14) 

(15) 


q,  =  ~ 


_d_  dT_  _  ST 
dt  Sx1  ^xi 


=  Lxxx  +  LxXx  + 


pAi&i 


pAx*i 

=  LiXx  +  pAx&x^a  -  2 
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and  the  upward  force  on  the  water  due  to  the  motion  of  the  case  is 


i  St  8t 

Qa  '  3F  aj '  '  ~i~ 


Note  that  the  force  is  not  simply  equal  to  an  added  mass  times  an  acceler¬ 
ation.  Velocity-dependent  terms  are  present  because  Lx  varies  with  xx 
and  xa. 

We  also  ignore  sound  waves  in  the  air  and  assume  that  the  air  is 
compressed  uniformly  and  isentropically ,  so  that  if  P  and  V  are  the 
instantaneous  pressure  and  volume  of  the  air 


PVY  =  PiV^  -  phvhy  -  povoy 


where  Px  and  Vx  are  the  initial  values,  P  and  V  are  the  values  that 
obtain  when  the  internal  pressure  is  equal  to  the  external  hydrostatic 
pressure,  Pq  and  Vq  are  the  values  that  obtain  when  the  piston  first 
touches  the  spring,  and  Y  =  1.4. 

The  "spring"  is  a  simple  energy-absorbing  device  whose  behavior 
approximates  an  ideal  linear  plastic  material.  By  ideally  linear  we  mean 
that  if  W(y)  is  the  work  done  in  dynamically  compressing  the  spring 
through  a  distance  y  (where  £  =  k3  -  £x),  then  the  resisting  force  is 


_  dW  . 
FS  *  37  *  ky 


where  the  stiffness  k  =  dsW/dy3  is  a  constant  independent  of  y  and  inde¬ 
pendent  of  the  rate  y.  By  ideally  plastic  we  mean  that  this  work  is  done 
irreversibly  so  that  if  the  compressive  force  were  relaxed,  the  compres- 
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sion  y  would  not  change.  In  practice,  the  "spring"  is  neither  perfectly 

linear  nor  perfectly  plastic,  but  these  imperfections  have  minor  effects 

which  can  be  analyzed  as  perturbations.  We  pre-position  the  spring  so 

that  contact  with  the  piston  first  occurs  when  V  ~  V„  and  P  ~  P„,  and 

o  H  o  H 

pick  a  spring  whose  stiffness  k  is  much  greater  than  the  equivalent 
stiffness  of  the  air  at  this  point.  In  this  way  the  impact  energy  is 
near  a  maximum,  the  stopping  distance  is  very  much  shorter  than  the 
travel  distance  (Vx  -  Vo)/A,  and  the  spring  force  is  the  dominant  force 
in  the  system. 

With  these  approximations,  the  equations  of  motion  for  the  piston 
and  case,  respectively,  are 

Ux*l  =  (P  “  Ph)Aj,  +  Fs  -  Qj.  (19) 

Ma52  =  -  (P  -  Ph)Ax  -  Fs  -  Qa  (20) 

where  Mx  and  Mg  are  the  masses  of  piston  and  case,  respectively.  The 
first  term  on  the  right  is  the  force  due  to  the  pressure  difference 
between  the  air  inside  and  the  water  outside;  the  second  term  is  the 
resistance  of  the  spring  and  is  present  only  when  the  piston  is  touching 
the  spring  and  compressing  it,  i.e,,  only  when  V  <  VQ  and  V  <  0;  the 
last  terms  are  the  reactions  on  the  piston  and  on  the  case  of  the 
unsteady  forces  in  the  water.  If  we  substitute  for  Qx  and  Qs  the  values 
given  by  Equations  (14)  and  (16),  we  have  the  equations  of  motion  in  a 
more  explicit  form 
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d  pAi^i 

^  [(Mx  +  Lx)^]=  (P  -  Ph)Ax  +  Fs  -  -y -  (21) 

,  pAx^x 

dt  [(Ms  +  L a)*3]  =  -  (P  -  Ph)Ax  -  Fg  +  -y-  (22) 

From  these  we  can  obtain  an  expression  for  the  conservation  of  linear 
momentum 


[(Mi  +  l>i)xx  +  (Mb  +  La)xa]  *  0  (23) 

Or,  since  &x  =  =  0  at  t  =  0,  then  for  all  t 

(Mx  +  Lx)ftx  +  (Me  +  Ls)&a  “  0 

and  this  relation  is  valid  whatever  the  form  of  the  internal  forces. 

There  is  also  a  geometrical  constraint 

$  -  Ax(*x  -  *3)  .  (25) 

Equations  (11),  (12),  (17),  (18),  (21),  (22),  and  (25)  are  suf¬ 
ficient  to  completely  determine  the  motions  Xx(t)  and  xa(t)  for  any 

prescribed  set  of  parameters  and  any  initial  conditions.  The  equations 
can  always  be  solved  on  a  computer  by  conventional  numerical  techniques 

to  any  desired  accuracy.  Furthermore,  Xj(t)  and  xa(t)  together  with 
Equations  (7)  or  (8)  completely  determine  the  magnitudes  of  the  radia¬ 
tion  multipoles  and  thus  completely  determine  the  sound  field. 

Thus  the  monopole  source  strength  is 
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Mi  +  Li 
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(26) 


where  the  reduced  mass  M  is 


I 

M 


1 

Mi  +  Li 


+ 


1 

M2  +  L2 


(27) 


Also  the  volume  acceleration  is 


••  Mi  +  Li 

V  = 


M 


Ai*xi  + 


Li&iAi 

Mg  +  L3 


■  ^  c*.  - 


A  *2 

pAiXx 
Mg  +  La 


.0 


Mi  +  Li  ^ 
*  M  ""  AlXl 


(28) 


regardless  of  the  nature  of  the  internal  forces.  The  peak  values  of  the 

necond  term  are  ordinarily  negligible  compared  to  the  first  term. 

To  calculate  the  dipole  radiation  we  take  the  origin  to  be  on  the 

longitudinal  axis  of  the  piston,  at  a  distance  bx  below  the  wetted 

surface  of  the  piston  and  a  distance  b8  above  the  bottom  of  the  case,  and 

take  n  as  a  unit  vector  pointing  upward.  Then  the  magnitude  of  the 
— o 

equivalent  dipole  force  is 

F  -  JJ  p^Crn*,)  d0  "  JJ  da  (29) 


■  p^lbxAi  +  pxa  ba As  -  Qi  -  Qa 
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r  /  Ml  +  LA  .  .  f  Ml  +  LAI 

-  x,  l^x  -  pA,bB  {  *  Lt  U,^ 

•  l1  *  ETii  -  sM 

-  «>(M,  +  10  [  CmT^TT  +  *7*u)  k‘  '  CmT^ZT  +  HTTu)] 


a  pA;  (Mg  -f  lip  ■  pApb^)  (M^  ’f  Li) 


M(M*  +  la) 


where  LQ  •  pAa(bl  +  ba)  depends  on  the  fixed  distance  bx  +  b,  and  is 
independent  of  the  position  of  the  origin  bt. 

Thus  the  amplitude  of  the  dipole  radiation  depends  on  bx  which 
defines  the  position  of  the  origin,  or  reference  "centroid,"  If  the 
added  inertia  Lx  were  independent  of  the  motion  (as  La  is),  i.e,,  if  tx 
were  zero,  and  if  the  centroid  is  taken  at  the  point  where 

Li  Lq  -  La 

Mi  +  Lj,  Ma  +  La 


+  Lj  Mj  *f  Lg 


then  the  dipole  force  F  would  be  identically  sero  (over  the  duration  of 
the  sound  pulse)  and  there  would  be  no  far-field  dipole  ratiation  relative 
to  this  centroid.  Since,  in  fact,  Lx  is  not  identically  zero,  there  is 
some  dipole  radiation  centered  at  this  centroid.  But  the  amplitude  of  the 
equivalent  dipole  nound  pressure  source  is  of  the  order  of  f  ■  77(pAiXi) 
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which  la  ordinarily  completely  negligible  compared  to  the  amplitude  of  the 
monopole  radiation.  Note  that  If  M»  •  M» ,  lM  »  Ltl  A\  »  A# ,  then  b*  -  b, 
which  meana  that  the  entroid  of  Kqnatton  (11)  la  halfway  between  the  top 
of  the  piaton  and  the  bottom  of  the  caae,  Alao  1 f  M»  ■*  U  »  N*  +  L* .  aa 
actually  occurs,  then  b^  *»  L^/nA)  and  is  independent  of  b», 

In  a  similar  manner  we  tan  show  that  the  quadrupole  radiation  ia 
negligible  for  this  simplified  implosion  mechanism  because  the  external 
and  internal  forcea  which  generate  the  surface  accelerations  have  no  net 
component  with  quadrupole  symmetry, 

We  can  astlmate  the  principal  features  of  the  sound  radiation 
quite  simply  if  we  neglect  the  force  (P  •  f^)A\  in  Equations  (19)  and 
(20)  in  comparison  with  tha  spring  forot  F  ,  Than  the  equation  of  motion 
during  Impact  simplifies  to 


1$  -f  ky  * 


(32) 


where  y  is  the  compression  of  the  spring  of  Equation  (18)  and  M  is  the 
reduced  mass.  The  maximum  value  of  the  term  on  the  right  ia  negligible 
compared  to  the  maximum  value  of  the  spring  force  ky  and  will  be  Ignored, 
The  initial  conditions  are  that  y  »  0,  J  ■  tuym  where  «*  «  k/M  end  is 
the  maximum  compreaaion  calculated  below  in  Equation  (36),  Hence  the 
solution  for  y  is  a  quarter  aine  wave 

V  ■  aPy  sin  tut  t  *  «/2to 

m 

y  »•  0  t  >  «/2cn  (33) 

and  the  monopole  sound  radiation  is  simply 
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Mi  +  l 


M 
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where  the  origin  ts  located  a  distance  bv  below  the  piston,  as  given  by 
Kquatton  (11), 

We  most  expect  that  eh*  necessary  deviation*  of  any  aotual  etrue- 
lure  from  thin  highly  simplified  modal  will  have  aoma  appreciable  effaces. 
These  wtll  be  discussed  when  we  deaertbr  the  actual  implosion  structure. 

We  also  derive  a  simple  estimate  of  the  maximum  compression  of 
the  plastic  spring.  Assume  that  at  the  time  of  maximum  compression  all 
parts  of  the  ayatem  are  tnatantaneously  at  rest  and  there  have  been  no 
energy  losses.  Then  W(yffl),  the  work  done  in  compressing  the  spring,  must 
equal  tho  work  done  by  the  hydrostatic  pressure  less  the  change  in  the 
internal  energy  of  the  air,  t,e, , 


(35) 


Thus,  if  W(v)  Is  known  independently  and  whatever  its  form  may  be,  then 

y  la  determined  by  the  Initial  parameters  P„  and  V  .  If,  in  addition, 
m  Ho 

we  specify  that  the  spring  is  linear,  then 


v  » .  2piv>r!ii  ( 

i  .  la v)  -  _L_ 

( Vi  V‘+  1 . 1 

"max  k  LP1  v 

Vv  /  +  V  ■  1  J 

O 

Note  that  these  relations  are  independent  of  the  masses  and  the  added 
masses. 
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RXPRRIMRNTAI,  VRRIPXCATXON  OP 
THR  THROl  I'lCAL  ANALYSIS 

RXPRR1MRNTAL  DRSXON 

Th«  Imploalon  device  uaed  In  the  experlmenta  waa  dealgned  Co  fit 
the  almpllfled  theory  and  to  ba  uaabla  at  deptha  up  to  1000  ft  (300  pal 
ovar  praaauraX  Tha  conatruotion  la  ahovn  in  Plguraa  2  and  3,  The  platon 
and  caaa  vara  tnada  of  atainlaaa  ataal  and  weighed  about  62  lb  in  air. 


Figure  2  -  Imploilon  Device  Reedy  (or  Tee* 
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Figure  3  -  Schematic  Drawing  of  Imploiion  Device 


The  upper  chamber  In  which  the  piaton  moved  hud  a  bore  diameter  of 
4 , 0 1  in,,  *  length  of  b,7  In,,  and  a  wall  thlckneaa  of  0,51  In.  Th#  upper 
chamber  was  Initially  cloaed  at  the  top  by  th#  platon  with  an  0-rtng  aenl. 
The  lower  chamber  contained  preamplifier*  for  th*  tnatrumenta, 

The  "breakbolt"  which  held  the  piaton  in  it*  Initial  position  waa 
screwed  to  the  piaton  at  the  bottom  and  pa*atd  through  a  hole  in  the 
narrow  "bridge"  Above,  This  braakbolt,  shown  in  Figure  4,  waa  made  from 
a  stock  high-strength  bolt,  1/2  in.  In  diameter  and  6  in.  long.  It  was 
necked  down  over  a  1/2- in,  length  to  a  diameter  which  waa  calculated  to 
break  suddenly  when  th*  outside  water  pressure  reached  some  preselected 
value.  Preliminary  testa  of  these  breakbolts  showed  that  the  effective 
tensile  strength  was  173  kpsl  ±  5  percent  and  that  the  elongation  before 
break  waa  leas  than  0,04  in. 

The  piaton  shown  in  Figure  5  had  a  diameter  of  4,00  in,  and  waa 
1.2  in,  thick  at  the  edge.  The  thickness  must  be  larger  than  the  maximum 
elongation  of  the  breakbolt  plus  the  contact  width  of  the  0-ring.  Also  a 
greater  thickness  minimises  the  possible  tilt  in  the  piston  as  it  travels 
down  the  Cube,  and  also  increases  the  frequency  of  the  lowest  elastic 
resonance  mode.  The  piston  was  honeycombed  by  drill  holes  in  order  to 
decrease  the  mass  without  a  proportional  decrease  in  the  stiffness.  Some 
auxiliary  experiments  showed  that  Che  lowest  elastic  resonance  mode  of  the 
piston  had  a  frequency  of  about  11,9  kH*  and  an  umbrella-like  vibration 
pattern  with  a  nodal  circle  whose  diameter  was  2,59  in. 

The  impact  "spring"  was  actually  a  ring  of  aluminum  (shown  in 

Figure  4)  which  was  annealed  dead  soft  to  make  it  plastic;  it  was  given  a 
triangular  cross  section  to  make  it  have  an  approximately  linear  response, 

i,e,,  to  make  the  force  resisting  compression  approximately  proportional 
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Figure  4  -  Braakbolti  and  Aluminum  Ring* 


PSD-330133 


Figure  5  -  Piston  with  Accelerometer  and  Balance  Weight 

to  the  amount  of  compression.  The  design  of  this  plastic  spring  was  sug¬ 
gested  by  experience  with  crusher  gages  using  soft  copper  balls. The 
diameter  of  the  impact  circle  of  the  aluminum  for  the  early  tests  was  3.0 
in.  (In  the  remaining  tests,  however,  it  will  be  made  equal  to  the  diame¬ 
ter  of  the  nodal  circle  of  the  vibration  pattern  of  the  piston  so  that 
impact  will  not  excite  the  fundamental  vibration  mode  of  the  piston.) 


Tlic  croM  itv cion  of  thw  aluminum  ring  was  in  the  shape  of  an  isosceles 
triangle  with  a  banc*  of  0,50  in,  and  a  height  of  0.50  in.  Thia  shape  waa 
chosen  bvcauao  auxiliary  teata  (described  in  Appandix  A)  ahowed  that  it 
had  a  reasonably  linear  response  with  an  acceptable  equivalent  stiffness 
of  k  »  3.98  x  I0fi  lb/in.  (for  the  3,0-in.  diameter).  The  aluminum  ring 
waa  cemented  to  a  pedeatal  cylinder  whoae  length  waa  selected  to  place 
cite  impact  at  any  dealred  point  along  the  chamber. 

An  accelerometer  waa  fastened  to  the  top  of  the  piston  at  a  point 
above  che  nodal  vibration  circle  in  order  to  measure  the  gross  vertical 
motion  of  the  piston  and  yet  be  insensitive  to  any  motion  in  the  funda¬ 
mental  elastic  mode,  (The  remaining  teata  will  use  two  identical  acceler¬ 
ometers  electrically  connected  in  parallel  and  placed  180  deg  apart  on  the 
nodal  circle  in  order  to  better  average  the  motion  over  che  face  of  the 
piston.)  The  accelerometer  was  a  piezoelectric  quartz-crystal  device 
(Kistler  Type  805A)  with  a  nominal  capacity  of  60  pF,  a  charge  sensitivity 
of  0.283  pC  per  g,  a  linear  output  range  up  to  108  g,  and  a  nominal  flat 
frequency  response  to  above  30  kHz,  The  accelerometer  was  connected  by 
2.5  ft  of  lo-nolse  coaxial  cable  (|i  Dot  50-3084)  to  a  preamplifier  (Ithaco 
Type  150  M102)  with  a  nominal  input  impedance  of  109  0,  an  output  imped¬ 
ance  of  50  U,  and  a  nominal  gain  of  unity.  A  small  capacitor,  which 
brought  the  total  input  capacitance  to  286  pF  was  placed  across  the  pre¬ 
amplifier  input  in  order  to  attenuate  the  expected  signal  below  the  peak 
linear  range  (3,8  V)  of  the  preamplifier.  The  accelerometer  case  was 
attached  to  the  piston  by  cementing  it  through  a  1-mil-thick  fiberglass 
base  with  Eastman  910  epoxy  adhesive;  it  was  further  insulated  from  the 
piston  and  the  water  by  a  coating  of  PRC  waterproofing  compound.  This  was 
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necessary  because  the  accelerometer  output  lead  is  electrically  grounded 
to  its  case. 


The  hydrophones  used  to  measure  the  sound  pressure  (Atlantic 
Research  Type  LC32)  are  made  of  three  cylinders  of  lead  zirconate-titanate 
each  1/4  in,  long  by  1/2  in,  diameter.  They  have  a  mechanical  resonance 
frequency  of  about  70  kHz,  a  capacity  of  about  15,000  pF,  and  a  nominal 
flat  frequency  response  to  about  20  kHz.  They  were  connected  by  3  ft  of 
lo-noiae  cabin  to  a  preamplifier  similar  to  that  used  for  the  accelerome¬ 
ter  (except  that  it  had  a  gain  of  five).  They  were  waterproofed  by  mold¬ 
ing  hydrophone  and  preamplifier  in  polyurethane  as  shown  in  Figure  6,  The 
sensitivities  of  the  hydrophone  and  preamplifier  combinations  were  deter¬ 
mined  as  described  in  Appendix  B;  they  ranged  from  10  to  40  pV/ubar. 

A  block  diagram  of  the  entire  instrumentation  is  shown  in  Figure  7. 
The  output  leads  of  the  preamplifiers  were  each  fed  through  a  long  cable 
(typically  Navy  Type  TTRS-4)  to  a  central  junction  box.  The  cable  also 
had  one  pair  of  leads  for  the  d-c  power  to  the  preamplifiers  and  another 


Figure  6  —  Hydrophone  ond  Preamplifier 
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Figure  7  -  Block  Diagrom  of  Instrumentation 

pair  of  leads  for  applying  a  voltage  calibration  to  the  input  circuit  of 
each  preamplifier.  The  signal  from  each  preamplifier  was  then  amplified 
by  an  Ithaco  Model  252  amplifier  whose  input  resistance  had  been  shunted 
to  about  50  fl  in  order  to  match  the  characteristic  impedance  of  the  long 
cable.  The  signals  were  then  recorded  with  a  tape  recorder  (Ampex  Model 
FR-1300)  operating  in  an  FM  mode  with  a  tape  speed  of  60  ips  and  an  upper 

frequency  response  limit  of  about  20  kHz,  A  precision  1-kHz  timing  signal, 

with  the  amplitude  coded  by  a  time  code  generator,  was  simultaneously 
recorded  on  one  channel  of  the  tape  recorder. 

After  the  signals  are  recorded,  they  can  be  visually  examined  in 
two  different  ways.  Individual  transient  pulses  can  be  examined  by  play¬ 
ing  back  one  channel  of  the  tape  through  a  CRO  and  triggering  the  sweep 
of  the  CRO  by  the  amplitude  of  the  transient  signal.  Or  all  channels  of 
a  tape  record  can  be  played  back  simultaneously  through  a  multichannel 
galvanometer-oscillograph,  e.g.,  Midwestern  Model  801  Visicorder,  running 
at  paper  speeds  up  to  64  ips.  The  available  galvonometers  have  a  flat 

frequency  response  to  about  5  kHz.  However,  if  the  magnetic  tape  is 
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played  back  at  7.5  ips  instead  of  at  the  recording  speed  of  60  ips,  the 
galvanometers  can  faithfully  record  input  signals  to  the  tape  with  fre¬ 
quency  components  greater  than  the  20-kHz  limit  of  the  tape  recorder.  A 
single  common  ground  for  the  entire  network  was  made  at  the  tape  recorder. 
The  net  frequency  response  limits  of  each  channel  are  discussed  in  Appen¬ 
dix  B. 

IMPLOSION  TEST 

The  experiments  were  originally  designed  to  be  conducted  over  a 
range  of  initial  conditions  in  deep  water.  However,  difficulties  in 
scheduling  the  facility  were  such  that  the  only  experimental  results  to 
date  (August  1970)  have  been  obtained  in  the  tank  of  the  Naval  Ordnance 
Laboratory  Hydroballistics  Facility. ^  This  100-ft-long  by  35-ft-wide 
tank  is  100  ft  deep,  but  at  the  time  of  the  test  it  was  filled  to  a  depth 
of  60  ft.  The  imploder  was  set  to  release  at  a  depth  of  40  ft,  and  three 
hydrophones  and  the  imploder  were  positioned  as  shown  in  Figure  8.  This 
arrangement  was  selected  so  that  (1)  a  sound  pulse  or  reasonable  ampli¬ 
tude  and  duration  would  reach  each  hydrophone,  (2)  these  pulses  would  not 
be  obscured  by  echoes  from  the  boundaries  of  the  tank,  and  (3)  comparable 
measurements  could  be  made  on  the  axis  of  the  imploder  and  in  a  direction 
of  90  deg  to  this  axis. 

The  three  hydrophones  were  each  suspended  at  their  nominal  posi¬ 
tions  by  1/8-in.  wire  rope,  with  the  hydrophones  pointing  in  a  vertical 
direction.  The  imploder  was  slowly  lowered  from  another  1/8-in  wire  rope 
at  a  rate  of  about  10  ft/min  until  the  breakbolt  fractured.  The  tape 
recorder  was  started  about  1  min  before  the  expected  fracture.  Voltage 
calibrations  were  imposed  on  the  preamplifiers  before  and  after  the  test. 
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After  an  implosion,  the  imploder  had  to  be  raised  in  order  to  replace  and 
position  the  breakbolt  and  the  aluminum  ring. 

For  the  test  at  a  40-ft  depth,  the  breakbolt  diameter  was  0.040  in. 
and  the  pedestal  height  was  4.02  in.  This  pedestal  height  wa3  estimated 
to  place  the  aluminum  ring  at  a  position  where  the  kinetic  energy  of  the 


HI,  H2.  H3  -  HYDROPHONES 
IMP  -  IMPLODER 


Pigur*  8  -  Location  of  imploder  and  Hydrophones  in  Tank  Test 
piston  and  water  would  be  a  maximum.  However,  it  is  not  critical  that 
these  pre-set  parameters  can  be  chosen  accurately  because  the  actual 
motion  of  the  piston  during  impact  is  measured  by  the  accelerometer. 

The  other  parameters  which  are  necessary  for  the  analysis  are 
listed  below.  The  asterisks  indicate  quantities  that  were  calculated 
for  the  actual  impact  position. 
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k  =  4.0  x  106  lb/in 

Mi  =  1.25  kg 

M2  =26.3  kg 

Li  =  0.90  kg 

L2  =  3.0  kg 

M  =  2.00  kg 

bj  =  3.0  in. 

=  44.3  ft 

P  =  2,30  atm 
H 

Pj  =  1.0  atm 

Vx  =  65.3  in3 

P  =  2.83  atm 
o 

V  =31.1  in3 
o 

In  this  test,  Lx  and  La  were  uncertain  by  perhaps  20  and  40  percent, 
respectively,  because  the  appropriate  coefficients  in  the  theoretical 
equations  (Equations  (11)  and  (12))  were  likewise  uncertain.  However, 
this  caused  a  smaller  uncertainty,  about  10  percent,  in  the  reduced  mass  M 
which  affects  the  amplitudes  of  the  acceleration  and  sound  pressure 
pulses,  an  uncertainty  of  about  5  percent  in  which  affects  the 
duration  of  the  acceleration  and  sound  pressure  pulses,  and  a  negligible 
uncertainty  in  the  ratio  (Mi  +  L^/M  in  Equation  (28). 

TEST  RESULTS 

The  teat  results  which  are  available  after  an  implosion  are:  (1)  a 
record  of  the  acceleration  of  the  piston  versus  time,  (2)  records  of  the 
sound  pressure  versus  time  at  each  hydrophone,  and  (3)  a  measurement  of 
the  set  deformation  of  the  aluminum  ring. 


Spring  stiffness 
Mass  of  piston 
Mass  of  case 

* 

Added  mass  for  piston,  Equation  (11) 
Added  mass  for  case.  Equation  (12) 

* 

Reduced  mass,  Equation  (27) 

* 

Centroid  position,  Equation  (31) 

* 

Depth  of  centroid  at  _mpact 

Hydrostatic  pressure 

Initial  air  pressure 

Initial  air  volume 

Air  pressure  at  impact,  Eq.  (17) 

Air  volume  at  impact 


25 


The  measured  compression  of  the  aluminum  was  .12, 6  ±  1,0  mils  which 
is  about  22  percent  below  the  value  of  41,6  mils  calculated  from  Equation 
(36),  The  difference  is  probably  due  to  energy  losses  as  the  piston  moves 
down  the  tube  and  to  the  fact  that  at  the  time  or  maximum  compression  some 
of  the  energy  must  have  gone  into  elastic  waves  in  the  structure.  In  fact, 
as  will  be  discuaaed,  the  accelerometer  record  can  he  Interpreted  to  ahow 
that  there  must  have  been  an  elastic  vibration  of  ths  piston  with  an 
amplitude  of  about  10  percent  of  the  plastic  displacsmsnt,  It  must  also 
be  verified  that  the  assumed  spring  stiffness,  which  was  measured  in  call* 
bration  tests  with  high- impact  energies  is  appllcsbls  for  ths  small  impact 
energy  of  this  tank  test. 

The  sequence  of  events  can  best  bs  verified  from  Figure  9  which  is 
a  Visicorder  record  of  the  1-kHs  timing  trace  and  ths  accelerometer  and 
hydrophone  signals.  The  first  signal  on  Che  accelerometer  record  wai  due 
to  the  impact  with  Lhe  aluminum  ring.  Prior  to  thia,  there  waa  no  aignal 
above  the  noise  level  at  the  expected  time  of  fracture  of  the  breakbolt, 
nor  any  signal  which  might  be  due  to  impact  of  the  moving  piaton  with  the 
side  of  the  case.  About  14  msec  after  the  impact  signal,  there  was  a 
sharp  pulse  of  unknown  origin;  this  may  have  been  due  to  mechanical  shock 
of  the  preamplifier  in  the  case.  The  signal  record  from  each  hydrophone 
shows  the  initial  sound  pulse  direct  from  the  implosion  followed  by  a 
series  of  echoes  from  the  top  (T) ,  bottom  (B) ,  or  sides  (S)  of  the  tank. 

Thus  Hi,  H2,  and  H3  identify  the  sound  pulses  direct  from  the  Implosion 
to  Hydrophones  1,  2,  and  3,  respectively;  Bl,  Tl,  and  Si  identify  the 
echo  pulses  which  reach  Hi  after  reflection  from  the  bottom,  top,  or  sides, 
respectively,  etc. 
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The  tabulation  given  below  indicate*  the  path  length  (aa  calculated 
from  the  geometry  of  Figure  d)  and  the  travel  time*  (aa  meaaured  in 
Figure  9)  for  each  sound  pulae  identified  in  Figure  9.  We  expected  that 
reflection  from  the  top  aurface  would  reverae  the  polarity  of  the  pulae, 
whereas  reflections  from  the  bottom  or  aidea  would  maintain  polarity. 
However,  the  apparent  echoea  from  the  bottom  and  aide*  do  nut  alwaya  ahow 
the  expected  polarity  (e,g,,  (83)  pulae*)  because  a  sheet  of  material  was 
hanging  vertically  in  the  tank  in  front  of  the  aide  wall,  and  this  material 
could  have  acted  aa  a  soft  reflector,  Some  pulaea  could  not  be  identified 
at  all.  Figure  10  la  a  plot  of  path  length  versus  travel  time.  The  data 
fit  a  straight  line  through  the  origin  with  a  slope  of  4830  fpa;  thle  la 
an  acceptable  value  for  the  speed  of  sound  in  city  water  at  a  temperature 
of  about  70  F, 
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We  ahall  anaiyaa  the  direct  pulaea  further;  however,  too  little  ta 
known  about  tho  nature  of  the  reflection  coefficient  from  each  aurface  to 
maka  it  profitable  to  examine  tha  aohoaa.  In  fact,  the  aehoaa  might  ba 
analysed  to  eatimate  tha  magnitude*  of  tha  reflection  coafficienta, 

Figure  U  ahowa  tha  accelerometer  aignal  with  an  acceleration  acalo 
datarminad  from  tho  voleaga  calibration  and  a  tima  acala  datarminad  from 
tha  timing  trace,  Tha  daahad  Una  curva  ia  tha  theoretical  curva 


«  10  40  M  M 

MTH  (.(NOTH  IN  PUT 

Figure  10  —  Path  Length  vertui  Travel  Tima  for  Hydrophone  Pulaa* 


Flgura  11  -  Platen  Acceleration  veraui  Tima 
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Vj  «  ^  L|'  «t“ym  Hitt  mt  I  m"  »  k/Ml  t  *  n/2/» 


with  y  taken  on  Ui  measured  value  (32,6  milt),  This  curve  has  been 
drawn  to  the  same  aoalaa  aa  the  measurtd  acceleration  and  la  arbitrarily 
placed  ao  that  the  peak  aocelarationa  occur  at  tba  earn#  time,  The  two 
curve*  are  groaaly  similar,  However,  the  peak  maaaured  acceleration  waa 
about  6  percent  above  the  theoretical  peak,  the  duration  of  the  maaaured 
pulae  waa  about  0,34  maec  compared  with  0,26  maec  for  the  theoretical 
curve,  and  the  meaaured  acceleration  had  an  oaciilatlon  at  about  IS  to 
20  kHv  auperimpoaed  on  the  mean  pulae.  These  differencea  could  all  be  due 
to  two  offsets  which  have  bean  ignored  in  the  simplified  theory:  (1)  the 
pis  ton  might  be  tilted  as  it  moves  down  the  tube  and  hit  the  aluminum 
spring  nonuniformly  and  (2)  there  are  undoubtedly  elastic  modes  of  the 
piston  and  case  with  periods  which  are  not  very  small  compared  to  the 
duration  of  the  pulae,  and  we  cannot  Ignore  them  by  analysing  the  piston 
and  case  as  rigid  bodies. 

Since  the  diameter  of  the  cylinder  bore  was  about  30  mils  more  than 
the  diameter  of  the  piston,  Chen  at  the  time  of  first  impact  with  the 
spring  the  center  of  the  piston  could  be  above  the  spring  by  as  much  as 


d  „  .vying _ x  30  -  hi  x  3f  -  40  mils 

edge  thickness  of  piston  1.1 


The  impact  must  then  rotate  and  straighten  the  piston,  as  well  as  deceler¬ 
ate  it.  The  duration  of  the  acceleration  pulse  will  be  increased,  but 
pending  a  more  detailed  analysis,  it  is  not  certain  what  effect  this  will 
have  on  the  amplitude  or  shape  of  the  acceleration  curve  at  a  single 


accelerometer.  Perhaps  the  fact  that  the  measured  compression  of  the  alu- 
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minum  was  fairly  uniform  around  the  rim  la  aoma  indication  chat  th#  p«ak 
deformation  and  paak  accalaratlon  are  not  greatly  affected,  and  perhapa 
chla  explain#  why  Che  peak  valuea  of  the  two  curvoa  in  Figure  11  agree  aa 
well  aa  they  do. 

There  are  aeveral  elaatic  mode#  which  could  explain  the  oacilla- 
tiona  on  the  accelerometer  record.  There  muet  be  a  longitudinal  vibration 
of  tho  aprlng  pedeatal  with  a  frequency  in  the  range  15  to  25  kHa,  There 
la  a  flexural  mode  of  Che  piaton  with  a  frequency  about  12  kHa.  There 
could  eaaily  be  elaatic  modea  in  the  caae  with  frequenciea  in  this  range. 
The  comparatively  alow  drop  In  the  acceleration  from  its  peak  value  and 
the  negative  peak  of  about  10  percent  of  the  positive  peak  must  be  due  to 
some  kind  of  elaatic  motion. 

We  next  compare  the  signals  from  the  three  hydrophones.  For  this 
purpose  they  are  shown  superimposed  in  Figure  12.  The  time  scale  has 
been  shifted  so  that  the  peaks  approximately  coincide  in  time,  and  the 
pressure  scale  has  been  magnified  in  direct  proportion  to  the  travel  dis- 
tance  from  implosion  to  hydrophone.  The  numerical  value  on  the  vertical 
scale  can  be  interpreted  as  the  sound  pressure  which  would  be  measured  at 
the  same  angular  direction  but  at  a  common  range  of  100  ft. 

There  was  good  agreement,  on  the  order  of  15  percent,  between  the 
three  records  for  the  first  0,3  msec  which  covered  the  development  and 
decay  of  the  principal  pulse.  After  0,3  msec,  the  signals  were  smaller 
and  differed  appreciably;  the  Hi  hydrophone  (on  the  longitudinal  axis) 
showed  the  most  complex  signal  and  the  H3  gage  (at  90  deg  to  the  axis) 
showed  the  simplest  signal.  This  makes  it  plausible  that  a  monopole 
source  dominated  in  the  first  0,3-msec  period,  whereas  a  longitudinal 
dipole  or  longitudinal  quadrupole  source  dominated  at  later  times.  This 


31 


Flgurt  12  —  Comparison  of  tht  Throo  Mydrophono  Rocoids  Adju»t«d  to  a  Common  Rang* 

is  in  accord  with  the  previous  discussion  of  the  acceleration  pulse 
because  (1)  the  rigid  body  rotation  of  the  piston  can  generate  no  monopole 
radiation  and  negligible  higher  multipole  radiation,  (2)  the  principal 
source  of  monopole  sound  is  due  to  the  plastic  motion  and  thus  should  be 
complete  within  0.3  msec,  and  (3)  the  longitudinal  elastic  motions  that 
persist  will  radiate  as  longitudinal  dipoles  or  quadrupoles. 

Since  the  sound  pressure  off  the  side  should  be  free  of  any  longi¬ 
tudinal  dipole  or  longitudinal  quadrupole  component,  we  compare  in  Fig¬ 
ure  13  the  sound  pressure  at  Hydrophone  3  with  a  theoretical  monopole 
radiation  at  this  range.  The  latter  was  calculated  from  the  volume 
acceleration  of  Equation  (28)  with  ^(t)  taken  as  the  experimental  values 
of  Figure  11,  Note  that  there  was  good  agreement  in  the  durations  of  the 
pulses  but  that  the  H3  signal  was  sharper  and  about  20  percent  higher  than 
the  calculated  monopole  radiation.  The  deviations  would  clearly  be  reduc¬ 
ed  if  we  could  correct  the  accelerometer  record  to  measure  the  mean  accel¬ 
eration  over  the  face  of  the  piston  rather  than  the  acceleration  at  one 
arbitrary  point. 
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The  experimental  data  are  too  aparae  to  justify  any  quantitative 
deduction  of  the  magnitudes  of  the  dipole  and  quadrupole  radiation  and  to 
explain  their  dependence  on  elastic  modes  in  the  imploder.  Likewise,  the 


Figure  13  -  Comparison  of  Sound  Pressure  Measured  at  Off-Side  Hydrophone  with  Theoretical  Me  topole  Component 

data  are  not  sufficient  to  justify  any  deduction  of  the  mean  motion  of  the 
piston  or  a  more  precise  calculation  of  the  monopole  radiation.  It  should 
also  be  added  that  the  hydrophone  sensitivities  are  not  known  to  better 
than  10  to  15  percent,  particularly  in  the  frequency  range  which  is  impor¬ 
tant  in  these  transient  pulses.  In  fact  there  is  some  small  additional 
evidence  that  the  hydrophones  may  appreciably  distort  sound  pulses  simi¬ 
lar  to  those  in  Figure  i2 ,  but  the  data  are  inconclusive. 

SUMMARY  AND  CONCLUSIONS 

1.  The  general  case  of  sound  production  by  an  imploding  structure 
is  analyzed  into  multipole  components,  where  the  amplitude  of  each  com¬ 
ponent  depends  on  a  particular  average  of  the  acceleration  and  pressure 
over  the  wetted  surface  of  the  structure.  In  most  cases,  the  dominant 
contribution  will  be  due  to  the  monopole  component  which  is  proportional 
to  the  second  time  derivative  of  the  volume  enclosed  by  the  wetted 
surface. 


33 


2.  In  the  particular  implosion  mechanism  used  in  the  experiments, 
the  principal  sound  production  is  associated  with  the  impact  of  two  parts 
of  the  structure  rather  than  with  the  initial  collapse  mechanism  or  with 
the  inrush  of  water. 

3.  The  theoretical  analysis  of  the  impact  mechanism  and  of  the 
generation  of  monopole  sound  radiation  is  substantially  verified  by  the 
limited  experimental  data. 

4.  There  are  significant  variations  in  the  experimentally  observed 
motions  and  sound  pressures  from  the  simplified  theoretical  analysis.  It 
is  plausible  that  these  are  due  to  nonsymmetric  impacts  and  to  the 
presence  of  low-frequency  elastic  vibrations  in  the  structure,  but  the 
experimental  data  are  too  sparse  to  establish  any  quantitative  relation. 

5.  For  the  additional  tests  planned,  the  implosion  depth  will  range 
up  to  500  ft,  and  the  sound  prer.**ure  will  be  measured  at  about  the  same 
range  at  points  above,  below,  and  off  to  the  side  of  the  implosion.  There 
will  be  more  accurate  measurements  of  the  surface-average  motion  of  the 
piston  and  more  accurate  voltage  calibrations  of  each  recording  channel. 
The  magnitude  of  the  elastic  oscillations  relative  to  the  plastic  motions 
should  be  reduced  because  the  support  pedestal  for  the  aluminum  will  be 
shortened,  the  piston  will  be  stiffened,  and  the  plastic  spring  will  be 
softer,  thus  lengthening  the  plastic  response  times  versus  the  elastic 
response  times.  It  is  expected  that  the  new  data  will  resolve  the 
uncertainties  in  the  present  analysis. 
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APPENDIX  A 


CALIBRATION  OF  THE  ALUMINUM  RINGS 

The  effective  dynamic  stiffness  of  the  aluminum  ring  was  determined 
by  calibration  tests  which  simulated  the  behavior  of  the  ring  in  the 
implosion. 

The  ring  was  placed  on  a  300-lb  steel  mass  which  rested  freely  on  a 
concrete  floor.  A  62.4-lb  brass  weight  was  dropped  from  an  adjustable 
height  and  impacted  the  aluminum.  The  falling  mass  was  guided  along  a 
0. 5-in. -diameter  rod  on  a  nylon  bushing.  Six  drop  tests  were  made  with 
two  different  sizes  of  rings.  One  ring  size  had  a  triangular  section 
1/2  in.  wide  by  1/2  in.  high,  and  the  diameter  of  the  impact  circle  was 
3,0  in.  The  size  of  the  other  ring  was  the  same  except  for  a  1,5-in.- 
diameter  impact  circle. 

Figure  14  shows  the  measured  values  for  the  mean  compression  of 
each  ring  after  impact,  versus  the  height  of  drop  divided  by  the  diameter 
of  the  impact  circle.  The  latter  coordinate  is  proportional  to  the 
impact  energy  per  unit  length  of  impact,  and  we  can  expect  this  parameter 
to  determine  the  resultant  deformation  so  long  as  the  initial  cross- 
section  shape  of  the  ring  was  the  same  in  ail  rings.  It  is  seen  that  all 
points  had  a  reasonable  fit  to  the  solid  line  with  a  slope  of  1/2,  thus 
indicating  that  the  final  deformation  was  proportional  to  the  square  root 
of  the  impact  energy.  It  is  true  that  the  2ta  fit  the  dashed  line,  with 
a  slope  of  0.51,  equally  well  or  better.  However,  the  square  root  rela¬ 
tion  is  within  the  experimental  error  and  is  easier  to  use  in  the 
analysis. 
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HEIGHT  OF  FALL/DIAMETER 

Figure  14  -  Calibration  of  Aluminum  Rings 


On  the  basis  of  this  solid  line  curve  and  assuming  that  the  bottom 
mass  does  not  move,  the  apparent  stiffness  of  the  3.0- in. -diameter  ring 
is 


k  -  (4.80  ±  .08)  x  10B  lb/ in. 
o 


(Al) 


However,  we  can  correct  this  apparent  value  for  the  effects  due  to 
the  finite  mass  and  the  probable  motion  of  the  bottom.  Thus,  let  and 
M2  be  the  top  and  bottom  masses,  respectively;  let  ux  be  the  impact  veloc¬ 
ity  of  the  top  mass  and  u2  be  the  common  velocity  of  the  top  and  bottom 
when  the  deformation  of  the  aluminum  reaches  its  maximum  value.  Then 
from  the  conservation  of  energy  and  assuming  that  the  bottom  is  free  to 
move 


Miui  Mi  +  Ma  k  a 

-  -  -  u3  + 

2  2 


(A2) 


and  from  the  conservation  of  momentum 
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Mlul  =  (M1  +  Ms)  u2 


(A3) 


Hence,  we  can  eliminate  us  and  find 


k  = 


Mg 

Mx  +  Mg 


(  =  (  ^  ^ 
V  y3  J  VMj  +  M2> 


k 

o 


(A4) 


and  for  the  aluminum  ring  with  a  3.0-in.  diameter 

k  =  3§§2~  x  4.80  x  10B  =  3.98  x  105  lb/ in.  (A5) 

The  proper  value  for  the  effective  stiffness  of  the  aluminum  ring  is 
somewhere  between  the  values  given  by  Equations  (Al)  and  (A5) ,  but  for 
the  high  frequency  motions  considered  here,  we  believe  that  the  founda¬ 
tion  was  effectively  soft  and  that  Equation  (A5)  is  the  more  accurate 
value. 


blank 
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APPENDIX  H 


SENSITIVITIES  AND  FREQUENCY  RESPONSE  CHARACTERISTICS 
OP  TRANSDUCERS  AND  ELECTRONIC  CHANNELS 

The  sensitivity  of  the  sccslsromstsr  was  checked  only  by  eompsrieon 
with  a  Keetler  Modal  S08K  accelerometer  which  la  alao  a  quartt-crystel 
device,  For  thie  purpoee  the  two  accelerometer*  vara  mounted  on  a  ahaka 
table,  with  tha  80SA  on  top  of  tha  808K,  Tha  frequency  waa  variad  over 
tha  ranga  30  to  1000  Ha,  and  tha  aignala  from  the  two  accelerometera  were 
compared  in  a  ■ imp la  tranafar  circuit,  Tha  aignal  amplitudes  ware  fourd 
to  ba  in  tha  ratio  of  tha  normal  aonaitivitiaa  within  about  2  percent. 
Hence  the  aanaitivity  of  the  SOSA  accelerometer  waa  taken  aa  0,283  pC 
par  g.  Note  that  both  accelerometera  are  aimple  quart*  crystal  devices 
which  can  be  expected  to  have  stable  and  predictable  responses  except  for 
accidental  mechanical  damage  or  accidental  electrical  leakage  to  ground. 

For  these  tests  it  was  desirable  to  have  an  absolute  calibration 
of  the  sensitivity  of  each  hydrophone  for  pulses  of  the  shape  shown  in 
Figure  12,  which  means  for  frequencies  in  the  range  0,1  to  20  kHz.  The 
calibrations  supplied  by  the  manufacturer  showed  a  variation  in  sensitiv¬ 
ity  of  the  order  of  ±  1,5  dB  over  this  frequency  range.  Furthermore,  it 
is  likely  that  the  sensitivity  at  any  one  frequency  in  this  range  is 
uncertain  by  1  or  2  dB.  Some  attempt  was  made  to  obtain  an  absolute  cali¬ 
bration  of  those  hydrophones  by  using  a  vibrating  column  of  water, ^  But 
this  method  could  be  used  only  at  low  frequencies  (up  to  400  Hz)  and  even 
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then  there  In  no  reliable  Indication  of  the  magnitude  of  the  absolute 
orror.  Accordingly,  the  aanaltivitlea  of  the  throe  hydrophones  were 
taken  aa 

Hydrophone  lj  7,6  pV/pbar 
Hydrophone  2 1  37,6  pV/pbsr 
Hydrophone  3:  38,0  pV/pbar 

Theae  are  the  manufacturer* '  value*  for  a  frequency  of  about  2  kH«,  We 
expect  that  theae  value*  may  be  off  by  a*  much  as  20  percent  In  absolute 
value.  However,  the  relative  sensitivities  of  two  hydrophones  are  prob¬ 
ably  accurate  to  perhaps  10  percent. 

A  voltage  calibration  of  any  accelerometer  or  hydrophone  channel 
is  best  made  by  imposing  a  known  step  voltage  (or  square  wave)  to  the 
calibration  resistor  in  the  input  circuit  to  the  preamplifier,  with  the 
transducer  connected  to  the  input,  and  with  the  amplifier  and  recorder 
gains  set  exactly  the  same  as  when  recording  the  signal.  Then  if  Rq  is 
the  calibration  resistor,  R  is  the  calibration  resistance  plus  any  lead 
resistance  plus  any  series  resistance  in  the  calibration  circuit,  C  is 
the  gage  capacitance  plus  any  lead  capacitance  plus  any  terminal  capaci¬ 
tance,  and  V  is  the  applied  step  voltage  across  R;  the  deflection  versus 
time  of  the  recorded  calibration  signal  is  exactly  the  same  as  though  a 
step  acceleration  &c(t)  were  imposed  on  the  accelerometer,  where 


a 

c 


SR 


(Bl) 
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and  where  S  ia  the  charge  aenaitivity  of  the  accelerometer,  Furthermore , 
cha  riae  time  and  ovarahoot  of  the  racordad  voltaga  atop  ia  a  good  meaaure 
of  tha  uppar  frequency  raaponaa  llmita  of  tha  antira  ampllflar  and 
racordar  ayatem,  whereaa  tha  dacay  In  tha  racordad  voltaga  atap  ia  a  good 
maaaura  of  tha  low* fraquancy  raaponaa, 

Figura  13  ahowa  tha  output  racord  of  a  10-Ha  aquara- topped  wave 
applied  to  tha  Input  of  tha  accalarometar  preamplifier,  The  output 
decayed  about  23  percent  In  30  maec,  Thla  lmpliea  a  decay  time  of  0.19 
aec  which  la  in  agreement  with  an  RC  of  the  Input  circuit  equal  to 
0,67  x  10° Q  timea  286  pF,  The  initial  part  of  the  rise,  shown  on  a  time- 
expanded  scale,  had  an  overshoot  of  about  15  percent  and  a  small  oscilla¬ 
tion  at  about  30  kHs.  Altogether,  this  voltage  calibration  showed  that 
a  voltage  pulse  will  be  recorded  with  negligible  error  provided  the  rise 
time  is  more  than  say  20  usee,  and  the  decay  time  is  less  than  1  msec. 

The  records  of  Figures  11  and  12  certainly  satisfy  these  criteria. 
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